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Abstract. The measurements of CMB anisotropy have opened up a window 
for probing the global topology of the universe on length scales comparable 
to and beyond the Hubble radius. For compact topologies, the two main 
effects on the CMB are: (1) the breaking of statistical isotropy in characteristic 
patterns determined by the photon geodesic structure of the manifold and (2) 
an infrared cutoff in the power spectrum of perturbations imposed by the finite 
spatial extent. We present a completely general scheme using the regularized 
method of images for calculating CMB anisotropy in models with nontrivial 
topology, and apply it to the computationally challenging compact hyperbolic 
topologies. This new technique eliminates the need for the difficult task of spatial 
eigenmode decomposition on these spaces. We estimate a Bayesian probability 
for a selection of models by confronting the theoretical pixel-pixel temperature 
correlation function with the COBE-DMR data. Our results demonstrate that 
strong constraints on compactness arise: if the universe is small compared to 
the 'horizon' size, correlations appear in the maps that are irreconcilable with the 
observations. If the universe is of comparable size, the likelihood function is very 
dependent upon orientation of the manifold wrt the sky. While most orientations 
may be strongly ruled out, it sometimes happens that for a specific orientation the 
predicted correlation patterns are preferred over the conventional infinite models. 



1. Introduction 

The remarkable degree of isotropy of the cosmic microwave background (CMB) 
points to homogeneous and isotropic Freidmann-Robertson- Walker (FEW) models 
for the universe. The underlying Einstein's equations of gravitation are purely local, 
completely unaffected by the global topological structure of space-time. In fact, in the 
absence of spatially inhomogeneous perturbations, a FRW model predicts an isotropic 
CMB regardless of the global topology. 

The observed large scale structure in the universe implies spatially inhomogeneous 
primordial perturbations exist which gave rise to the observed anisotropy of the CMB. 
The global topology of the universe does affect the local observable properties of the 
CMB anisotropy. In compact universe models, the finite spatial size usually precludes 
the existence of primordial fluctuations with wavelengths above a characteristic scale 
related to the size of the universe. As a result, the power in the CMB anisotropy is 
suppressed on large angular scales. Another consequence is the breaking of statistical 
isotropy in characteristic patterns determined by the photon geodesic structure of the 
manifold. One can search for such patterns statistically in the COBE maps, and to 
the extent that they are not there, one can constrain the size of the universe and its 
topology. 
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Much recent astrophysical data suggest the cosmological matter density 
parameter, f2 , is subcritical jy. If a (possibly varying) cosmological constant is 
absent or insufficient to bring the total density to the critical value, this would 
imply a hyperbolic spatial geometry for the universe (commonly referred to as the 
'open' universe in the cosmological literature). The topologically trivial (simply 
connected) hyperbolic 3-space, Ti. 3 , is non-compact and has infinite size. There are 
numerous theoretical motivations, however, to favor a spatially compact universe Q. 
To reconcile a compact universe with a flat or hyperbolic geometry, consideration of 
spaces with non trivial topology (non simply connected spaces) is required. A compact 
cosmological model can be constructed by identifying points on the standard infinite 
flat or hyperbolic FRW spaces by the action of a suitable discrete subgroup, T, of the 
full isometry group, G, of the FRW space. The infinite FRW spatial hypersurface is the 
universal cover, tiled by copies of the compact space (most appropriately represented 
as the Dirichlet domain with the observer at its basepoint) . Any point x of the compact 
space has an image = 7^x in each copy of the Dirichlet domain on the universal 
cover, where 7^ G Y. 

The hyperbolic manifold, H 3 , can be viewed as a hyperbolic section embedded 
in four dimensional flat Lorentzian space. The isometry group of Ti 3 is the group of 
rotations in the four space - the proper Lorentz group, 5*0(3, 1). A compact hyperbolic 
(CH) manifold is then completely described by a discrete subgroup, T, of the proper 
Lorentz group, 5*0(3, 1). The Geometry Centre at the University of Minnesota has 
a large census of CH manifolds and public domain software SnapPea ||. We have 
adapted this software to tile H 3 under a given topology using a set of generators of 
r. The tiling routine uses the generator product method and ensures that all distinct 
tiles within a specified tiling radius are obtained. 

A CH manifold, M, is characterized by a dimensionless number, Vm = Ym/oI 3 , 
where Vm is the volume of the space and d c is the curvature radius [[| . There are a 
countably infinite number of CH manifolds with no upper bound on Vm ■ The smallest 
CH manifold discovered so far has Vm = 0.94 g. [| The Minnesota census lists several 
thousands of these manifolds with Vm up to ~ 7. In the cosmological context, the 
physical size of the curvature radius d c is determined by the density parameter and the 
Hubble constant H : d c = (c/Hq)/t/1 — Sl ■ The physical volume of the CH manifold 
with a given topology, i.e., a, fixed value of VM/d 3 , is smaller for smaller values of Slo ■ 
Two quantities which characterize linear dimensions of the Dirichlet domain are i?> 
and i?<, the radii of circumscribing and inscribing spheres, respectively. 

In the standard picture, the CMB that we observe is a Planckian distribution 
of relic photons which decoupled from matter at a redshift IIOO. These photons 
have freely propagated over a distance R LS ~ 2d c arctanh-y/1 — fl , comparable to 
the "horizon" size. For the adiabatic fluctuations we consider here, the dominant 

contribution to the anisotropy in the CMB temperature measured with wide-angle 

1/2 

beams (9 FWHM comes from the cosmological metric perturbations through 

the Sachs- Wolfe effect. 

The adiabatic cosmological metric perturbations can be expressed in terms of a 
scalar gravitational potential, $(x, r). The dynamical equation for the gravitational 

X The volume of CFL_manifolds is bounded from below and the present theoretical lower bound 
stands at V^v) > 0.167 pi. There exist sharper lower bounds within subclasses of CH manifolds under 
restrictions on topological invariants |?|. It has been conjectured that the smallest known manifold 
is in fact the smallest possible |5|. 
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f{t) = ! ^3 ' . (1) 



potential allows for separation of the spatial and temporal dependence, | ^(x, r) = 
i 7, (r)$(x, t ls ), where F(t) encodes the time dependence of the metric perturbations 
and $(x, r LS ) is the field configuration on the three- hypersurface of constant time 
r = t ls when the last scattering of CMB photons took place. We shall study open, 
CIq < 1, models with zero cosmological constant where in the matter dominated 
phase, § 

5(sinh 2 t — 3r sinh r + 4 cosh r — 4) 
(coshr — 1) 

Here and further on we use dimensionless conformal time r expressed in units of the 
curvature radius. A non-zero cosmological constant can be trivially incorporated in 
our analysis by using the appropriate solution for F(t). 

We write the Sachs- Wolfe formula for the CMB temperature fluctuation, AT(q), 
in a direction <j, in the form 

AT 1 /" x h d 

— (q) = -Hq Xll ,r LS ) + 2 d X f( X )HqX,r w ) , /(x) = -^(r) , (2) 

where x is the affine parameter along the photon path from % = at the observer 
position to x H — R hs /d c . The first term is called the surface or "naive" Sachs- Wolfe 
effect (NSW). The second term, which is nonzero only if $ varies with time between r LS 
and now, is the integrated Sachs- Wolfe effect (ISW). The angular correlation between 
the CMB temperature fluctuations in two directions in the sky is then given by 

C(q,q') = /^L(q)^L(q>)\ = I , r LS )$(g'x H , r LS )> 



dx fix) [(^(qx,T LS )^{q'x H ,T Lfi )) + (Hq'x,n s )${qxv,n s ))] 

H d Xl fix,) I H dx 2 /(%) (^(q Xl ,T LS )^(q'x 2 ,T LS )) . (3) 
Jo 

The main point to be noted is that C(q, q') depends on the spatial two point correlation 
function, = ($(x, t ls )&(x' , t ls )) of <i> on the three-hypersurface of last scattering. 
This is due to the fact that the equation of motion for $ allows a separation of spatial 
and temporal dependence. 

Although in this work we restrict our attention to the Sachs- Wolfe effect which 
dominates when the beam size is large, we should point out that other effects which 
contribute to the CMB anisotropy at finer resolution can also be approximated in terms 
of spatial correlation of quantities defined on the hypersurface of last scattering B . 



2. Method of Images 

As described in the previous section, the angular correlation function for the CMB 
temperature anisotropy can be expressed in terms of the two point spatial correlation 
function £$(x, x') for the gravitational potential field, $(x, t ls ) defined on the three- 
hypersurface of last scattering. The correlation function can be expressed formally 
in terms of the eigenfunctions, of the Laplace operator, V 2 , on the hypersurface 
(with positive eigenvalues k\ > 0), as J8|, 

^(x,x / ) = ^*i(x)**(x')P$(/c i ), where (V 2 + fc 2 )* 4 =0. (4) 

i 

§ At the scales appropriate to CMB anisotropics, damping effects on $ can be neglected. 
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The function P$(fc^) is the initial power spectrum of the gravitational potential. 

In order to calculate the on CH manifolds we have developed the method of 
images | 1 1 which evades the difficult problem[| of solving for eigenfunctions of the 
Laplacian on these manifolds and, in general, on any nonsimply connected manifold, 
M.. The manifold M. can be regarded as a quotient of a simply connected manifold 
M u by a fixed point free discontinuous subgroup T of the isometry group of M u , 
i.e., M. = Ai u /T. The manifold M. u is referred to as the universal cover of A4. The 
elements 7 S T produce images of the manifold M which tessellate M u . 

The spectrum of the Laplacian on a compact space (thus with closed boundary 
conditions) is a discrete ordered set of eigenvalues {k 2 } (k 2 = and k 2 < kf +1 , with 
multiply repeated eigenvalues counted separately). The set {^i} is the corresponding- 
complete orthonormal set of eigenfunctions on A4. The eigenfunctions ^ on M. 
are also eigenfunctions on the universal cover, M. u . As a result these eigenfunctions 
simultaneously satisfy the integral equations for the correlation functions in both Ai u 
and M, 



f dx' a(x,x') ^(x') = PtfcMx). 
f dx' ££(x,x') *i(x') = P 9 (ki)Hfi(x) . 

J M 



(5) 
(6) 



Using the fact that the are automorphic with respect to T (\P(7x) = ^(x), V7 € T), 
and that M. tessellates M u , 

f dx' a(x,x , )* l (x') = y; / d* 1 a(x )7 xO 



7er' 



dx' 



M 



7er 



^(x')> (7) 



where denotes a possible need for regularization at the last step when the order 
of integration and summation is reversed. The regularization is achieved through 
a counterterm which removes the contribution of the homogeneous k 2 = mode. 
Using eqs. (0) and ((?]) leads to the main equation of our method of images Q 
expressing the correlation function on a compact space (and more generally, any non- 
simply connected space) as a sum over the correlation function on its universal cover 
calculated between x and the images 7X' (7 £ T) of x': 



a(x, x ') = £a(x,7x') 



7er 



^a(x,7x')-^/ dx'a(x,x') 



(8) 



Obtaining closed form expressions for eigenfunctions may not be possible beyond the sir 



the 



topologies. Examples where explicit eigenfunctions have 
and noncompact hyperbolic space with horn topology 

geodesic flow on CH manifolds presents an inherently difficult problem even in terms of a numerical 
estimation of eigenfunctions. 



beenused include flat models jl2,|l3| |14|. 
3, Il5| . The generically chaotic nature of 
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3. Compact Hyperbolic Models 

For cosmological CH models, M. u = Ti. 3 , the three dimensional hyperbolic (uniform 
negative curvature) manifold. The local isotropy and homogeneity of H? implies 
£|(x, x') depends only on the proper distance, r = d(x, x'), between the points 
x and x'. The eigenfunctions on the universal cover are of course well known for 
all homogeneous and isotropic models p7| . Consequently £$(x, x') can be obtained 
through equation (|J). The initial power spectrum P$(k) is believed to be dictated by 
an early universe scenario for the generation of primordial perturbations. We assume 
that the initial perturbations are generated by quantum vacuum fluctuations during 
inflation. This leads to 



where j3 = ^/{kd c ) 2 - 1 and P # (/3) = /3(f3 2 + l)P$(fc) / '(2tt 2 ) is the mean square 
fluctuation in <E> per logarithmic interval of k. We consider the simplest inflationary 
models, where P$(&) is approximately constant in the subcurvature sector ( kd c > 1). 
This is the generalization of the Harrison- Zeldovich spectrum in spatially flat models 
to hyperbolic spaces |l8[ . Sub-horizon vacuum fluctuations during inflation are not 
expected to generate supercurvature modes, those with kd c < 1, which is why they 
are not included in cq.(||). Indeed, since H 2 > l/(ad c ) 2 , for modes with kd c < 1 we 
always have kj (aH) < 1 so inflation by itself does not provide a causal mechanism for 
their excitation. Moreover, the lowest non-zero eigenvalue, ki > in compact spaces 
provides an infra-red cutoff in the spectrum which can be large enough in many CH 



spaces to exclude the supercurvature sector entirely {k\d c > 1). (See § [hl| and M.) 
Even if the space does support supercurvature modes, some physical mechanism needs 
to be invoked to excite them, e.g., as a byproduct of the creation of the compact space 
itself, but which could be accompanied by complex nonperturbative structure as well. 
To have quantitative predictions for P$(fc) would require addressing this possibility 
in a full quantum cosmological context. We note that our m ain c onclusions regarding 



peculiar correlation features in the CMB anisotropy (see § 4.2 ) would qualitatively 
hold even in the presence of supercurvature modes. 

3.1. Numerical implementation of method of images 

Although equation ^ encodes the basic formula for calculating the correlation 
function, it is not numerically implementable as is. Both the sum and the integral in 
equation (^|) are divergent and the difference needs to be taken as a limiting process 
of summation of images and integration up to a finite distance r* : 



£|(x,x') = lim 



rj <r. 



= d(x,7 i x'). (10) 



The volume element in the integral corresponds to Ti 3 . In || we conduct an instructive 
analytical study of the regularization procedure by varying an artificial infrared cutoff. 

In hyperbolic spaces, the number of images within a radius r» grows exponentially 
with r* and it is not numerically feasible to extend direct summation to large values 
of r*. The presence of the counterterm, however, besides regularizing, significantly 
improves convergence. This can be intuitively understood as follows : ^( r j) represents 
a sampling of a smooth function at discrete points rj. In a distant radial interval 
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[r, r+dr], r 3> i?>, dr ~ i?> there are approximately (Att/Vm) sinh 2 r images. The 
sum, J2 r ^*( r j')' w hhin this interval is similar to the (Monte-Carlo type) estimation 
of the integral, therefore one may approximate the sum over all distant images beyond 
a radius r* by an integral to obtain 

A />0O 

a(x,x') = £ a(r 3 ) + t- dr sinh 2 r aW . (H) 

r~^r, V M Jr, 

The tilde on £J(x, x') denotes the fact that it is approximate and unregularized. 
Subtracting the integral (Att/Vm) dr sinh 2 r £|(r) as dictated by the regularization 
eq.fll), we recover the hnite r* term of the limiting sequence in eq. (fTo)) . This 
demonstrates that even at a finite r*, in addition to the explicit sum over images 
with rj < r* , the expression for £J in eq. ( ^o|) contains the gross contribution from all 
distant images, r* > r*. Numerically we have found it suffices to evaluate the right 
hand side in eq. (|lQ) up to r* about 4 to 5 times the domain size i?> to obtain a 
convergent result for £J(x, x'). 

After this project was completed, we learned that estimation is similar to the 
remainder term obtained in |2C| to incorporate the gross contributions of long periodic 
orbits in evaluating the "Selberg trace formula" (e.g., |2l| ) in the quantum particle 
problem on 2D CH surfaces. Whereas the focus in |20) was an accurate computation 
of eigenvalues for determining the quantum energy levels of the classically chaotic 
system, our focus is on correlation function evaluation in 3D. We emphasize that it 
is the regularized procedure jl0| ) following directly from eq. (^) which is central to the 
success of our method. The form of the unregularized correlation function, £J(x, x'), 
in eq.([y]) suggested by the approximate treatment of distant images discussed above, 
in fact, also follows from the form of the regularizing counterterm in eq.(||). 

Figure [l] illustrates the steps involved in implementing the regularized method 
of images. The value of £|> as a function of r* has some residual jitter, which arises 
because of the boundary effects due to the sharp "top-hat" averaging over a spherical 
ball chosen for the counterterm in eq. ©. This can be smoothed out by resummation 
techniques |22|. We use Cesaro resummation for this purpose. 



3.2. Computing the CMB anisotropy 

For Gaussian perturbations, the angular correlation function, C(q 1 q') 1 completely 
encodes the CMB anisotropy predictions of a model. To confront observations, 
the celestial sphere is discretized into N p pixels labeled by p, with N p determined 
by the angular resolution of the measurements. We test the models against the 
four-year COBE-DMR data. The six COBE-DMR maps [^3| are first compressed into 
a (A+B) (31+53+90 GHz) weighted-sum map, with the customized Galactic cut 
advocated by the DMR team, basically at ±20° but with extra pixels removed in 
which contaminating Galactic emission is known to be high, and with the dipole and 
monopole removed. Although one can do analysis with the map's (2.6 ) 2 pixels, this 
"resolution 6" pixelization of the quadrilateralized sphere is oversampled relative to 
the COBE beam size, and there is no effective loss of information if we do further 
data compression by using "resolution 5" pixels, (5.2°) 2 (24|]. The celestial sphere is 
then represented by N p = 1536 pixels before the Galactic cut, with N p = 999 pixels 
remaining after the cut is made. 

Full Bayesian linear comparison of the model with the data requires computation 
of the Np x N p pixel-pixel correlation matrix, Ct P p' = C(q p ,q p i). The expression for 
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Figure 1. This illustrates the regularization of the correlation function £|,(x, x'). 
The example shown is the correlation function at zero-point separation at some 
point on the CH manifold m004(-5,l). The topmost panel shows the contributions 
to ££(x, x') from close images and the distances to the images. The upper 
solid curve in the middle panel shows the cumulative build up of < r ?*( r j ) 
with successive addition of distant images. The dashed curve is the regularizing 
counterterm required to remove the zero-mode contribution and the lower solid 
curve shows the cumulative value of regularized §^(x, x') which fluctuates around 
the true value after sufficient number of images have been added. The bottom 
panel demonstrates that this residual jitter (s°) in the estimate of £Jj,(x,x') can 
be removed by Cesaro resummation. The sequence of lines, s 1 ,s 2 & s 3 shows the 
result of applying first, second and third order Cesaro resummation. The accuracy 
at the second order is usually sufficient. 



C(<7> q') in eq. (Q) involves an integral along the line of sight from the observer to the 
surface of last scattering - the integrated Sachs- Wolfe term (ISW). We find that a 
simple integration rule using Nl points along each line of sight gives accurate results 
when the points are spaced at equal increments of F(t), which is determined by SIq. 
Consequently, to get the full Ct pp ' matrix we need to evaluate the correlation function, 
£a,(x, x'), between N p (N p + l)iV£/2 pairs of points on the constant-time hypersurface 
of last scattering. 

To calculate Ct pp > in the standard infinite open cosmological models with this 
real-space integration to an accuracy comparable to that of the traditional evaluation 
in k-space, Nl ~ 10 is sufficient. Moreover, the real space integration is much faster 
in terms of CPU time. For the method to remain accurate in compact models, Nl 
should, of course, exceed the number of the times a typical photon path crosses the 
compact Universe. We found Nl ~ 10 is still enough for the models we have analyzed 
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so far. 



4. Results 



In 8 4.1 



In 6 4.2 



we discuss the power spectrum of fluctuations in CH spaces, 
discuss how the tessellation of H 3 by the finite domains is reflected in the CMB 
anisotropy correlation function. In § 4.3, we present the results of full Bayesian 



probability analyses of large angle CMB anisotropy predictions for a selected set of 
CH models. 



4-1- Power spectrum 

When going to a compact space from its universal cover, the power which is spread 
out in a continuum in fc-space gets redistributed and accumulates in peaks at the 
discrete eigenvalues allowed in the compact space. For example, in the well known 
case of T 3 (flat 3-torus), the power spectrum is a series of delta functions at 
k = 2TTy/ (n x /L x ) 2 + (n y /Ly) 2 + (n z /L z ) 2 , where n x , n y and n z are integers. 

By contrast, for CH manifolds precise spectra of eigenvalues of the Laplacian 
are not known. We can obtain some information about the spectrum by applying 
the method of images to each (3 mode in the integrand for the correlation 
function at zero-lag £|,(x, x). The quantity V$(k, x) estimated in this way is 
defined by £|(x, x) = L d\n(k)V^,(k,x) and describes the contribution to £§,(x, x) 
from a unit logarithmic band dln(fc) = 1. For a discrete spectrum, V^(k.x) = 
J2i 5(k — ki)V<s,(k) \^j(k, x)| 2 (here, differently from eq.(Q), i enumerates only the 
distinct eigenvalues and j labels the degenerate eigenfunctions belonging to the same 
eigenvalue; S is the radial delta- function) . Indeed, in the formalism of images the 
singular delta-functions in the spectrum are recovered, in principle, by the precise 
cancellation of the contributions from all images. Determination of the spectrum 
V^{k), in this sense, is a far more difficult task than that of the correlation function 
which is non-singular. Our approximation technique includes gross integral estimation 
of the impact of distant images and results in a spread-out convolved spectral power 
distribution. Increasing r* progressively sharpens spectral profiles near true positions 
of discrete eigenvalues |25| . 

In general, CH spaces are not only globally anisotropic but also globally 
inhomogeneous, even though the universal covering space (here, 7i 3 ) is homogeneous 
and isotropic. For V^{k,x) in compact space the inhomogeneity manifests itself in 
the dependence on x. Figure |^ shows a sample of power spectra for some of the 
CH manifolds at two random positions x. There is a definite signature of a strong 
suppression of power at small /3 in all the cases that we have explored. Since this 
suppression is observed at all x, we conclude that the small j3 part of the spectrum 
is absent in the models considered. This is qualitatively similar to the infrared cutoff 
known for the compact manifolds with flat and spherical topology. Quantitatively too, 
the break appears around k ~ 0{RZ}), consistent with intuitive expectations. 

An infra-red cutoff at the lowest non-zero eigenvalue, ki > 0, exists for all compact 
spaces. The Cheeger's inequality p6fl provides a lower bound on k\ for a compact 
Riemannian manifold, M: 

2 hs~* ^.(iS") 
fc 1 > T , hc=M min{v[Mi)iV[M2) y (12) 
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where the infimum is taken over all possible surfaces, S, that partition the space, 
M, into two subspaces, M x and M 2 , i.e., M = M x U Af 2 and 5 = <9A/i = dM 2 (S 
is the boundary of Mi and M 2 ). The isoperimetric constant he depends more on 
the geometry than the topology of the space, with small values of he achieved for 
spaces having a "dumbbell-like" feature - a thin bottleneck which allows a partition 
of the space into two large volumes by a small-area surface 26, E7j. More regular 



compact spaces do not allow eigenvalues which are too small. For example, the Cheeger 
limit for all flat T 3 manifolds is k\ > 2/L, where L is the longest side of the torus. 
Although estimation of he is not simple, for any compact space, M, with curvature 
bounded from below there exists a lower bound on he in terms of the diameter of M , 
cLm = s ^P x ,yeM d>{x,y) (2^, |21J1 ; for a 3-dimensional CH space, 

-l 



fci > hc/2 > 



1 



I'M 



1/2 



dt cosh (t) 



= 0.92/d, 



(13) 



This result prohibits supercurvature modes for all CH spaces with 4m < 0.92 d c . 

There are also upper bounds on k\. The bound, k\ < Ahc/d c + I0h c , |^9| 
does not allow for a firm conclusion that the space supports supercurvature modes, 
kid c < 1, unless (using eq. (|l3|)) 4m > 10.6 d c (i?> > 5.3 d c ). Upper bounds based on 
comparison with the first Dirichlet eigenvalue on a subdomain of A4 cannot impose 
k\d c < 1, since Dirichlet eigenvalues cannot be less than d' 1 . 



4-2. CMB anisotropy correlation function 

In standard cosmological models based on a topologically trivial space, such as Tt 3 , 
the observed CMB photons have propagated along radial geodesies from a 2-sphere 
of radius i? LS (that we refer to as the sphere of last scattering, SLS) centered on 
the observer. The same picture also applies to CH models when the space is viewed 
as a tessellation of the universal cover. If the compact space fits well within the 
SLS, the photons propagate through a lattice of identical domains. As a consequence, 
strong correlations build up between CMB temperature fluctuations observed in widely 
separated directions. The correlation function C{q,q') is anisotropic and contains 
characteristic patterns determined by the photon geodesic structure of the compact 
manifold. These correlations persist even in CH models whose Dirchlet domains are 
comparable to or slightly bigger than the SLS. This is the key difference from the 
standard models where C(q,q') depends only on the angle between q and q' and 
generally falls off with angular separation. 

In Figure^, the complex behavior of the correlation function is illustrated with the 
example of the small CH model (m004(-5,l)). The SLS encompasses « 150 domains 
for flo = 0.3 and ~ 20 domains for Slo = 0.6 and is comparable to the size of one 
domain for f^o = 0.9. The angular correlation along any arbitrary great circle in the 
sky in lower f2o models shows distinct peaks as one encounters repeated copies of the 
Dirichlet domain. The peaks are more pronounced when one considers only the surface 
terms of eq. (||) in the CMB anisotropy. Including the line of sight ISW contribution 
tends to smear the peaks but adds its own characteristic features |^| . 

Consider compact universe models which are small enough such that the SLS does 
not completely fit inside one domain, R LS > R < . The CMB temperature is expected 
to be identical along pairs of circles if temperature fluctuations are dominated by the 
surface terms at the SLS |50|. We identify these matched circles on the sky in our 
models and check the extent of cross correlation seen at the angular resolution of 
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Figure 2. The primordial power spectrum Pj^(0) exhibits suppression of long- 
wave power if the Universe has compact topology. P£ (/3) = 1 corresponds to an 
infinite open Universe. Examples of the spectrum cutoff in several CH Universes 
are shown. The first 6 have about the same value of ij> / d c , where i?> is the 
minimal radius of the sphere which completely encompasses the Dirichlet domain. 
These spaces have volumes V jd\ around unity. The sixth space is m004(-5,f), and 
the seventh one is the v3543(2,3) space with volume Vm/^c = 6.45, which are 
used as examples of "small" and "large" spaces in this paper. The plots show the 
integrand for the correlation function £|,(x,x') at zero separation, x = x', at two 
different locations (left and right panels). This contrasts with isotropic spaces for 
which the integrands will be the same. The vertical line at k = 2/iJ> marks the 
value below which there is strong suppression of large-scale power . 



COBE-dmr. Figure || show the matched pairs of circles in two ( "small" and "large" ) 
CH models superimposed on random realizations of the theoretical sky generated 
from the computed pixel-pixel correlation matrix for these models. The models were 
chosen to have a volume comparable to the volume V SL s within the SLS. Even with the 
coarse pixelization of COBE-dmr, we do see fairly good cross-correlation in the CMB 
temperature along matched circles in our realizations. Again, the pattern of correlated 
circles is more pronounced when the surface terms are dominant, as happens in the 
small model (SCH) with f2 = 0.9. The cross-correlation coefficient p between the 
circles is in the range p as 0.6 — 0.95 in this case, whereas it is in the range p as 0.2 — 0.6 
for the large model (LCH) with ilo = 0.6. 

The circles of identified pixels on the SLS are not the whole story H . Enhanced 
NSW cross-correlations, but at a lower level, also exist between all pairs of points 
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Figure 3. This shows the behavior of the correlation function of the CMB 
temperature on a great circle in the sky in a CH model (m004(-5,l)). The solid 
curves in the first and third panels show C(6) for the surface term in the Sachs- 
Wolfe (NSW) effect for Cl = 0.3 k, 0.6, respectively. They reflect the spatial 
correlation along a circle on the sphere of last scattering (SLS). The number of 
peaks in C(0) matches the number of Dirichlet domains that the circle intersects. 
The dashed curves show corresponding results for the simply connected infinite H 3 
models with Qq = 0.3 &z 0.6. The second and fourth panels are analogous to the 
first and third panels, where the integrated Sachs- Wolfe (ISW) effect is included. 
In this more physically correct case the sharp NSW peaks have been diluted by 
the ISW contributions. However, ISW induces new features, in particular, note 
the appearance of strong negative correlations in panels two and four. 



on the SLS which are projected close to each other on the CH manifold. This 
can be seen as secondary maxima in the examples in Figure^; these are absent in 
standard cosmological models. Some features persist at a detectable level even when 
the compact universe encompasses the SLS, i.e.,R LS < i?< , although circles are 
absent in this case (the effect dies out for R LS <C R<). When the relative ISW 
component is significant, the geometrical patterns based on pointwise identifications on 
the SLS are supplanted by more complex features arising from identifications between 
photon geodesies, e.g., the strong negative correlations evident in Figure||. Correlation 
features also arise from the global inhomogeneity of CH spaces since it implies that 
the rms temperature fluctuations varies with location in the sky in a pattern set by 
the topology of the space. 

The reader's intuition can be sharpened with examples from simpler topologies. 
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Flat torus models exhibit forward-backward correlation symmetry along their axes J12| , 
|l4j . Patterns in CMB temperature maps have also been discussed for other flat models 
in jl5| . An example of a hyperbolic model is the non-compact horn-like space |l6| for 
which |l5| argue the CMB sky will exhibit a flat spot. 

4-3. Bayesian analysis constraints from COBE— DMR 

As discussed above, compact models which are not much larger than the sphere of last 
scattering tend to show significant correlated patterns in the CMB sky. To the extent 
that these correlations are absent in the COBE-DMR data, the models tend to be at 
odds with observations. Figures || and |^ compare theoretical realizations of the CMB 
anisotropy in the LCH and SCH models with the COBE-DMR data. They should be 
compared with the 'DATA' map in Fig. ||, a Wiener filtered picture of the CMB data, 
using a theoretical model which best-fits the data. What one should be noting is the 
shapes of the patterns and not the specific locations of the patterns, since these can 
change from realization to realization. The full Bayesian analysis takes into account 
all possible realizations. The incompatibility of models with small Vm/Vs LS (SCH- 
£lo = 0.3, 0.6) is visually obvious: the best fit amplitudes are high which is reflected 
in the steeper hot and cold features. Although, SCH-f2 = 0.9 and LCH-£1 = 0.6 do 
not appear grossly inconsistent, the intrinsic anisotropic correlation pattern is at odds 
with the data statistically. 

On each of our selection of models, we have carried out a fully Bayesian analysis 
of the probability of the model given the COBE-DMR 4yr data. In Table |l|, we present 
the relative likelihood of the selected models to that of the infinite H? model with the 
same SIq- (The COBE data alone does not strongly differentiate between the infinite 



Table 1. The Log-likelihoods of the compact hyperbolic models relative to the 
infinite models with same f2o are listed. The probabilities are calculated by 
confronting the models with COBE-DMR data. The values quoted are likelihoods 
marginalized over the amplitude of the initial power spectrum. The volume within 
the sphere of last scattering (SLS) relative to the volume of the compact models 
of the universe (Vsls/^m) is listed. The three columns of logarithm of likelihood 
ratios C/Co correspond to the best, next best and the worst values that we have 
obtained amongst 24 different rotations of the compact space relative to the sky. 
The number v in brackets gives the conventional, albeit crude, translation of the 
probabilities to a Gaussian likelihood C/Cq ~ exp[— u 2 /2]. Only the last model 
in a specific orientation appears to be consistent with the COBE-DMR data. 



CH Topology 


n 


Vsls/Vm 


Log of Likelihood Ratio (Gaussian approx.) 


[V M ,R>/dc] 








Orientation 










'best' 


'second best' 


'worst' 




0.3 


153.4 


-35.5 (8.4<t) 


-35.7 (8.4ct) 


-57.9 (10.8o-) 


m004(-5,l) 


0.6 


19.3 


-22.9 (6.8<x) 


-23.3 (6.8<r) 


-49.4 ( 9.9<t) 


[0.98,0.75] 


0.9 


1.2 


-4.4 (3.0cr) 


-8.5 (4.1<t) 


-37.4 ( 8.60-) 


v3543(2,3) 


0.6 


2.9 


-3.6 (2.7a) 


-5.6 (3.3<t) 


-31.0 (7.9a-) 


[6.45,1.33] 


0.8 


0.6 


2.5 (2.2ct) 


-0.8 (1.3<t) 


-12.6 (5.0o-) 



hyperbolic models with different Slo-) The anisotropy of the theoretical C(q, q') causes 
the likelihood of compact models to vary significantly with the orientation of the 
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space with respect to the sky, depending on how closely the features in the single 
data map available to humans match (or mismatch) the pattern predicted in C(q, q'). 
Some optimal orientations may also have the "ugly" correlation features hidden in the 
Galactic cut. We analyzed 24 different orientations for each of our models and found 
that only the model with Vm > Vsts (LCH-O = 0.8) cannot be excluded. (At one 
orientation this model is even preferable to standard CDM; this raises a question of 
the statistical significance of any detection of intrinsic anisotropy of a space when only 
a single realization of data is available.) 

Similar conclusions were reached by some of the authors (JRB, DP and 
I. Sokolov [l41 ) for flat toroidal models. Comparison of the full angular correlation 
(computed using the eigenfunction expansion, eq.(|J)) with the COBE-DMR data led 
to a much stronger limit on the compactness of the universe than limits from other 
methods jl2[ [l3|. The main result of the analysis was that V sls /Vm < 0.4 at 95% CL 
for the equal-sided 3-torus. For 3-tori with only one short dimension (or for the 
non-compact 1-torus), the constraint on the most compact dimension is not quite as 
strong because the features can be hidden in the "zone of avoidance" associated with 
the Galactic cut. 

Although our results strongly indicate that manifolds with small Vm/V sls are 
unlikely to survive confrontation with the COBE-DMR data, we emphasize that when 
this ratio becomes of order unity very interesting correlation patterns appear which, for 
certain orientations, may even be preferred by the COBE-DMR data. Since manifolds 
with large VmIo% allow low values of Slo while still avoiding small Vm/V sls , even large 
mean curvatures may be compatible, and certainly the f2o ~ 0.3 — 0.6 values that some 
of the astrophysical data currently prefers can still be nicely accommodated within 
the CH framework, ^j] Of course, there are many more manifolds for which R < may 
be small but i?> large, for whom the constraints will be even more dependent upon 
how specific orientations may deliver unwanted correlation features into the zone of 
avoidance, as the asymmetric 3-torus results reveal, or line up predicted correlation 
patterns with chance upward or downward fluctuations in the data. When R < is large 
compared to R hs , we expect the results will quickly converge towards the usual infinite 
hyperbolic manifold results. The intermediate terrain still encompasses ample scope 
for interesting topological signatures to be discovered within the CMB. Although our 
methods are quite general, testing all manifolds in the SnapPea census this way is 
rather daunting, and there are countably infinite manifolds not yet prescribed. What 
may be more promising for discovery are specialized statistical indicators, which are 
less sensitive than the full Bayesian approach we have used here, but not as manifold 
sensitive; e.g., the statistical techniques which exploit the high degree of correlations 
along circle pairs that |50| have emphasized, and Fig. |] reveals. Maps like we have 
constructed will be necessary to test the statistical significance of such methods. We 
also note that dramatically increasing the resolution beyond that of COBE-DMR is 
quite feasible with current computing power using our techniques. 

CH models with non-zero cosmological parameter, A, and associated d ensi ty. C2a are certainly 
possible. The COBE-DMR data by itself does not strongly constrain Oai e.g., jlOj, and the constraints 
presented here would roughly correspond if Qq is replaced by Qo + Q\, although there will be be 
quantitative differences. 
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Figure 4. This shows two full-sky CMB anisotropy maps, plotted as pairs 
of 180° diameter hemispherical caps, one centered on the South Galactic Pole 
(SGP) and one on the North (NGP). They are one of an infinite number of 
possible random realizations based on computed pixel-pixel correlation matrix for 
the model in question. Both surface and integrated Sachs- Wolfe effects have been 
included. The power was normalized to best match the COBE data. In contrast to 
Figs, and q the maps are not optimally filtered. The model labels L(arge)CH 
and S(mall)C H refer to the CH models v3543(2,3) and m004(-5,l), respectively. 
(The model number associated with the topology corresponds to that of the census 
of CH spaces from the Geometry center, Univ. of Minnesota; SCH is one of the 
smallest and LCH is one of the largest spaces in the census). The value of f2o in 
each was chosen so that Vm ~ V'sls ■ The matched pairs of circles expected if the 
CMB anisotropy is dominated by the surface terms |3o| are superimposed on the 
map for each model. Each pair is labeled by the same number centered on the 
circles (and the same colour). The relative phase is shown by identified points 
marked by a diamond and a triangle on each circle in a pair. For clarity, we show 
only the eight largest pairs out of 35 for the LCHq — q.q case. Even at COBE-DMR 
resolution, we find the cross-correlation between the temperature along matched 
circles is very good in the SCHQ 0= o.g model. The ISW contribution is larger 
at f2rj = 0-6, and the cross-correlation coefficients are systematically smaller for 
LC-Hq 0= o.6 circle pairs. The contours are linearly spaced at 30 fiK steps. 



Figure 5. The figure consists of a column of three CMB sky- maps showing a 
pair of 140° diameter hemispherical caps each, centered on the South (SGP) and 
North (NGP) Galactic Poles, respectively. The top map labeled DATA, shows the 
COBE-DMR 53+90+31 GHz A+B data after Wiener filtering assuming a standard 
CDM model, normalized to COBE. The next two maps are one random realization 
of the CMB anisotropy in v3543(2,3) - our choice of a L(arge)CH model example, 
for f!rj = 0.6 & 0.8 based on our theoretical calculations of C(q, q') convolved with 
the COBE-DMR beam. Both surface and integrated (ISW) Sachs- Wolfe effects have 
been included in C(q, q'). No noise was added. The power was normalized to best 
match the COBE data. The theoretical sky was optimally filtered using the COBE 
experimental noise to facilitate comparison with data. LCH with Qrj = 0.8 is 
compatible with the data with a suitable choice of orientation while f2rj = 0.6 is 
ruled out (See Table [l]). For all the maps in Figs. ^Jand|iJ the average, dipole and 
quadrupole of the \b\ > 20° sky were removed and a 20° Galactic latitude cut was 
used, with extra cuts to remove known regions of Galactic emission proposed by 
the COBE team accounting for the ragged edges. The contours are linearly spaced 
at 15 fiK (AT/T = 0.55 ■ 10" 5 ) steps. The maps have been smoothed by a 1.66° 
Gaussian filter. 



Figure 6. The three CMB sky-maps showing a pair of 140° diameter 
hemispherical caps each, centered on the South (SGP) and North (NGP) Galactic 
Poles, are analogous to the lower two plots in Fig. pmut for the S(mall)CH model 
m004(-5,l). The fact that f2o = 0.3 & 0.6 models above are strongly ruled out 
by the COBE data is obvious visually, while the f2o = 0.9 can be excluded on the 
basis of our Bayesian analysis. (See Table hi) 
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